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1. Introduction. The phenomenological shell model remains the bedrock of nuclear 
structure physics Il|,@]. In its standard form, we accept the empirically established notions 
of closed shell or magic nuclei and of single particle or single hole excitations with respect to 
these special cores (below, for brevity, we speak only of particles). To study the properties of 
nuclei which are removed from these stable cores by two or more nucleons, one adds residual 
two-particle forces. To understand low energy behavior of low and medium mass nuclei 
one restricts the allowed single-particle excitations and residual interactions to the valence 
shell. The solution of the resulting matrix diagonalization problem, which is straightforward 
in principle, has been achieved only up to mass number A^ = 48 |^ because of the rapid 
growth of the dimensionality of the Hamiltonian matrices. Beyond that there are several 
possibilities: One can study somewhat heavier nuclei with Monte Carlo calculations @] 
that utilize the entire valence shell-model space. One can also reduce the dimensionality 
of the Hamiltonian matrices to tractable sizes in several ways, either by utilizing only the 
lowest irreducible representations of relevant approximate symmetries 0,^ or by applying 
the variation- Al method to a trial space suggested by the deformed shell model 0,0] . Except 
for the first example cited 0, these latter methods allow one to break the bounds of the 
valence shell restriction. 

More than three decades ago, A. K. Kerman and one of the authors proposed an alterna- 
tive to the standard linear approach to the shell model. Originally designed as a method for 
restoring the broken symmetry of mean field solutions [Q , it soon became clear that it was 
a general formulation of quantum mechanics |T^ that could also be used to study the shell- 



model problem. It was argued that especially in cases of well-developed collective motion, 
one could replace the linear methods that use a large basis of states by a non-linear method 
involving a tractable set of states. Early attempts to apply this method to semi-magic nuclei 
n|,0 were at best only modestly successful and were not followed up. 



Our aim in this letter is to reactivate the original program by reporting a successful 
application to a deformed nucleus in the rare earth region, ^^^Gd. For such a nucleus 
not only is it technically impossible to apply the valence shell model based on spherical 
single-particle excitations, but the restriction to the valence shell itself fails badly ||l|,P-§. 
Starting from a spherical shell model expanded to include all orbits bound in a realistic 
(Wood-Saxon) single-particle potential and a standard Hamiltonian widely applied for heavy 
nuclei, we describe a fully microscopic derivation of some of the properties of the ground state 
rotational band including energies, charge and mass quadrupole matrix elements, and pairing 
matrix elements. This work was carried out using results obtained from a semi-microscopic 
description, described below and referred to as CPC, of the low energy properties of odd 
deformed nuclei. The results reported for ^^^Gd are almost certainly not special to this 
nucleus, nor is the method necessarily confined to the ground-state band. 

2. Model and method. As Hamiltonian we choose the form 



H = Y1 haittlitt. 
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In the first term a and a^ are the spherical shell model annihilation and creation operators, 
a labels the principal and angular momentum quantum numbers of states in a spherical 
Wood-Saxon potential with spin-orbit coupling, the subscript a the same set minus the 
magnetic quantum number, i distinguishes neutron from proton, and hai are the eigenvalues 
in the respective wells. The second term is a sum of products of mass multipole moments 
with strengths Kij-i] in the following we shall retain only the most important terms, those 
with L = 2, though in a more refined treatment, we should include L = 4 |T^. The last term 
is a sum of pairing interactions with strengths Qi-i , of which we include only the dominant 
monopole, though here the L = 2 term can also be considered to be well established [[7|. 
The adequacy of the Hamiltonian ([^) as a representation for those properties of a more 
realistic interaction that lead to collective behavior has been carefully documented in a 



recent investigation [^ 



To explore the consequences of (|TD, we calculate the commutator of a and of a^ with the 
Hamiltonian and take the matrix elements between states \J^iy) of a chosen odd nucleus 
of mass number A^, with angular momentum quantum numbers J, /i and all other labels 
indicated by u and the corresponding states of the relevant even neighbors, \IMn{N ± 
1 >. Suppressing charge quantum numbers, we encounter in the single-particle terms the 
coefficients of fractional parentage (CFP) 

Vj^^{a; IMn > = < JiJiv\a^\IMn{N + 1) > (2) 

Uj^,{a- IMn) = < Jfiu\al\IMn{N - 1) > . (3) 

For the evaluation of a typical interaction term, consider, in an abbreviated notation, 

(J|ag2|/(A^ + 1)) 

= Y.Vj{n{I'{N + l\Q,\I{N + l)), (4) 

/' 

which involves only the completeness relation. By this means the non-linear terms are 
expressed as sums of products of terms in which one factor, the CFP, depends on the odd 
nucleus, whereas the other depends on the properties of the even cores. With a corresponding 
treatment of the pairing interaction, we obtain equations with characteristics and properties 
that we now describe. 

In addition to the CFP, which are coupled by the pairing interactions, there occur in 
these equations matrix elements of the mass quadrupole moments for two different, neighbor- 
ing, even nuclei and matrix elements of the pairing interaction between the two neighbors. 
Together we refer to these as the core matrix elements. The structure of these equations. 



which is given in all detail in [|15[ and will not be reproduced here, bears a striking resem- 
blance to those of the Hartree-Bogoliubov mean-field theory, but in contrast to the latter 
our equations are formally exact, conserving both angular momentum and particle number. 
If we assume that the core matrix elements are known, the resulting equations are linear 
and define an Hermitian eigenvalue problem for the energies of the odd nucleus relative to 
the average ground ground state energy of the neighboring cores. In this interpretation, the 
chemical potential of the odd nucleus and the excitation spectra of the even neighbors are 
added to the list of quantities assumed to be known. The solutions are mutually orthog- 
onal, and the normalization will be considered below. The possibilities inherent in such a 



generalized semi-microscopic theory, first noted in ||T6[, was first developed and applied by 



Donau and Frauendorf [17,18|. Recently further development and applications have been 



carried out by the authors in a series of papers of which the latest are [p!3i p!5| . 

It is vital to recall the approximations that are made in order to make even the linear 
scheme workable. The most important can be understood by examination of @. If the 
starting state of the core, \In) belongs to a given low-lying rotational band, then we know 
that intra-band transitions are by far the dominant ones, though a few neighboring bands 
provide some residual strength, and these are included in the calculations in order to satisfy 
the sum rule. Once the choice of core bands has been made, there remains a vital question 
associated with the space of the single odd nucleon. For all the examples done, we find that 
results for the observables of interest have essentially converged when three major shells are 
included. For the purposes of the new results reported in this letter, we have nevertheless 
done calculations that include all bound orbitals. The reason for this will be explained 
below. 

3. Self consistency: particle number. We now take the next step and ask whether the 
solutions of the CPC calculations that we have described are in fact self consistent. Naturally 
we have chosen to examine first the most favorable case of a strongly deformed nucleus in 
a region where the properties to be tested vary slowly and smoothly with particle number. 
We start with what may at first sight appear to be a trivial example, the conservation of 
particles. The operator for the total particle number, 

can be separated into a sum of four terms 

N = iVp,+ + iVp,_ + iVn,+ + iVn,_ , (6) 

where the subscripts distinguish charge of the nucleons and parity of the single-particle 
orbitals. Because we include in the Hamiltonian (|I|) only multipoles of even parity, each 
of these quantities is conserved. The even nucleus chosen for study plays the role of the 



heavier of the two cores in a calculation ||T^ carried out for Gd, an axially deformed 



nucleus with states \IMK{N + 1)), where K is the angular momentum of the band head. 
Further discussion will be confined to the ground-state band with K = and this quantum 
number will be suppressed. We thus need the four sets of eigenvalues 

iV,,± = Y.(^M{N + l)|aL±a,,±|/M(iV + 1)) 

= E v};m-Jm)v}^mjm), (7) 

each of which should be independent of /, as it automatically is independent of M. To 
evaluate these sum rules, one needs the CFP for the neutron levels of ^^^Gd of both parities, 
which had been obtained earlier [jl9| and the proton levels of ^^^Eu, which were obtained for 
present use. The results are shown below in Table |. It has been verified that the sums in 
(^) depend on I only in the third decimal place. 

We consider the results given in the table to be strongly encouraging. In this regard, 
two points must be noted. The first is that we must finally confront the problem of the 
normalization of the CFP. Though we failed to emphasize the point in our previous work. 



transition matrix elements calculated in CPC are independent of an overall rescaling of 
the normalization, provided it is the same for all states \J^v). In our work we assumed 
unit normalization, as in the strong coupling core-particle model, but this choice cannot be 



exact, as we know from our early work on spherical nuclei [^. To make the appropriate 
corrections requires incorporating into our algorithm a set of sum rules derived from the 
Fermion anti-commutation relations. This has not yet been done. Second, to achieve a 
result so close to the exact one, we must include (numerically significant) contributions from 
a large number of solutions of the eigenvalue problem, including high-lying ones that play 
no role in the fit to the known observables of the odd nuclei. These points are relevant as 
well for the remainder of our discussion but will not be mentioned again. 

4. Self consistency: quadrupole matrix elements. In the reference Hamiltonian (|l]) there 
occur three quadrupole coupling constants. In the case under discussion, the experimental 
electric quadrupole matrix elements are in good agreement with rigid rotor values. For 
the CPC calculations we then assume proportionality between neutron and proton mass 
quadrupole elements, as expressed by the relation 

{IMK\Q''\I'M'K') = 7] {IMK\QP\I'M'K'), (8) 

where ?7 is a constant evaluated below. It follows that for the core-particle theory, we can 
work with the proton quadrupole moment alone provided we introduce different effective 
coupling strengths for the neutron and proton spectra, according to the equations 

l^f = (l^PP + V /^pn) , (9) 

Kf = (ri Atnn + /^pn) • (10) 

Thus in fitting CPC to the data in the odd nuclei, we are allowed to choose and indeed 
find slightly different values for the effective coupling constants. For these purposes and for 
the further development, the actual value of r] is reflected only in the value that has to be 
assigned to the effective coupling strengths. 

Since we are dealing with operators of the form 

their core matrix elements are again quadratic sums in the same set of CFP as enter the 
calculation of the number. The first test of self-consistency is that these sums have a "shape" 
consistent with the rigid rotor assumption, as expressed in (||), a test that is passed with 
flying colors. This is seen partly from Fig. |l|, which emphasizes the fact that not only is the 
angular momentum dependence of the electric quadrupole matrix elements given correctly, 
but also their magnitudes and signs. From the fact that the neutron quadrupole matrix 
elements follow parallel curves, we deduce the value r] = 1.1. 

5. Self consistency: pairing. For the pairing matrix elements, which are linear combina- 
tions of matrix elements of type 

{IMn{N - l)\aaaa\I'M'n{N + 1)) 

= E U},M; IMn)Vj,,{a- I'M'n'), (12) 

JflU 



there are again two tests. First there is the requirement that the matrix elements be inde- 
pendent of angular momentum, as was assumed in the input. Second we must reproduce 
the value of this matrix element. For the neutron pairing, the diagonal matrix elements vary 
monotonically between 1.859 for / = and 1.814 for 1 = 8, compared with the experimental 
value of 1.65 for the gap parameter. The corresponding values for the proton pairing are 
1.672 for / = and 1.612 for 1 = 8, compared with the experimental value of 1.321. 

6. Self consistency: moment of inertia. There remains the test of the self consistency 
of the excitation spectrum. This requires first that the diagonal elements of (^ satisfy the 
rigid rotor equation 

Ej = {IM\H\IM) = ^^^^^, (13) 

and second that the experimental value of X be reproduced. The first of these requirements 
is well-satisfied, and therefore we can confine our attention to the moment of inertia. In 
Table |T[ we display not only the final calculated value of the moment of inertia, but also 
the contributions of individual terms of H, broken down according to nucleonic charge. The 
absence of contributions from the neutron quadrupole moment simply reflects the fact that 
we organized the calculation so that the quadrupole term is expressed completely in terms the 
proton quadrupole operator and of the effective coupling strengths. Notice that the major 
contribution comes overwhelmingly from the single-particle term and that the contribution of 
the quadrupole term is insignificant. (This term then contributes only deformation energy.) 
The self consistency is as close as one has a right to expect, in view of the well-known effect 
of the quadrupole pairing interaction |^, which is not in our calculation. 

7 Final remarks. The method of choosing parameters used in fitting the odd-nucleus 
data has been thoroughly covered in our previous work |T^|T^JT^. In the present work 



we have checked a number of sum rules that should be satisfied by a reasonably complete 
CPC calculation and found agreement to within 10 percent. These results imply that for 
a limited number of states we have demonstrated a new route for passing directly from 
a spherical shell model to the properties of a deformed rare-earth nucleus. The immediate 
next steps are twofold: to add to the working Hamiltonian the additional simple interactions 
mentioned in the text and to add to the algorithm a proper formulation for normalization 
of the coefficients of fractional parentage. 
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TABLES 



(iV) 


Actual 


Calculated 


neutrons + 


44 


45.06 


neutrons — 


52 


50.64 


neutrons total 


94 


95.70 


protons + 


38 


36.26 


protons — 


26 


28.90 


protons total 


64 


65.16 



TABLE L 



Particle number, actual and calculated, for ^HCd. 



E2 — Eq 


Tsp 


Iq 


Xa 


2total 


-'-cxp 


neutron + 


0.00291 


- 


0.00121 






neutron — 


0.00532 


- 


-0.00236 


- 




neutron total 


0.00823 


- 


-0.00157 


- 




proton + 


0.001943 


0.00141 


-0.00232 


- 




proton — 


0.001272 


0.000523 


0.00194 


- 




proton total 


0.003215 


0.001933 


-0.00126 


- 




Total 


0.01145 


-0.000643 


0.001144 


0.011951 


0.0124 



TABLE IL Contributions to the moment of inertia arising from the different terms in the 
Hamiltonian. Here sp refers to the single particle contribution, Q to the quadrupole contribution, 
and A to the pairing term. 
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FIG. 1. Values of the matrix elements of Qp for ^^^Gd. The negative values represent the 
diagonal elements {IK = 0\\QP\\IK = 0) whereas the positive values the off-diagonal elements 

{IK = 0\\QP\\I + 2K = 0). 



